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The width of the giant dipole resonance (GDR) at finite temperature T in 120 Sn is calculated 
within the Phonon Damping Model including the neutron thermal pairing gap determined from the 
modified BCS theory. It is shown that the effect of thermal pairing causes a smaller GDR width 
at T < 2 MeV as compared to the one obtained neglecting pairing. This improves significantly the 
agreement between theory and experiment including the most recent data point at T = 1 MeV. 
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I. INTRODUCTION 



Intensive experimental studies of highly-excited nuclei during the last two decades have produced many data on 
the evolution of the giant dipole resonance (GDR) as a function of temperature T and spin. The data show that the 
GDR width increases sharply with increasing T from T > 1 MeV up to ~ 3 MeV. At higher T a width saturation has 
been reported (See for the most recent review). The increase of the GDR width with T is described reasonably 
well within the thermal shape-fluctuation model and the phonon damping model (PDM) [J- [3- The thermal 
shape-fluctuation model assumes an adiabatic coupling of GDR to quadrupole degrees of freedom with deformation 
parameters j3 and 7 induced by thermal fluctuations and high spins in the intrinsic frame of reference. Although this 
model shows an increase of the GDR width with T comparable with the experimental systematic at 1.2 < T < 3 
MeV, the GDR shapes generated using the strength function of this model differ significantly from the experimental 
ones Q . The PDM considers the coupling of the GDR to pp and hh configurations at T ^ as the mechanism of the 
width increase and saturation. The PDM calculates the GDR width and strength function directly in the laboratory 
frame without any need for an explicit inclusion of thermal fluctuation of shapes. The PDM reproduces fairly well 
both of the observed width HEl and shape 0,0 of the GDR at T ^ 0. 

In general, pairing was neglected in the calculations for hot GDR as it was believed that the gap vanishes at T = T c < 
1 MeV according to the temperature BCS theory. However, it has been shown in [j| that thermal fluctuations smear 
out the superfluid-normal phase transition in finite systems so that the pairing gap survives up to J 1 > 1 MeV. This 
has been confirmed microscopically in the recent modified Hartree-Fock-Bogoliubov (HFB) theory at finite T [jj, 
whose limit is the modified-BCS theory [HI El- Other approaches such as the static-path approximation shell 
model calculations 01 > as wei l as the exact solution of the pairing problem [l4| also show that pairing correlations 
do not abruptly disappear at T =/= 0. It was suggested in |l 51 that the decrease of the pairing gap with increasing 
T, which is also caused by pp and hh configurations at low T, may slow down the increase of the GDR width. By 
including a simplified T-dependent pairing gap in the cascade calculations using the PDM strength functions, Ref. 
has improved the agreement between the calculated GDR shapes and experimental ones. 

Very recently the 7 decays were measured in coincidence with 17 particles scattered inelastically from 120 Sn [l6j. 
A GDR width of around 4 MeV has been extracted at T = 1 MeV, which is smaller than the value of ~ 4.9 MeV for 
the GDR width at T = 0. This result and the existing systematic for the GDR width in 120 Sn up to T ~ 2 MeV are 
significantly lower than the prediction by the thermal shape-fluctuation model. Based on this, Ref. 0] concluded 
that the narrow width observed in 120 Sn at low T is not understood. The aim of the present work is to show that 
it is thermal pairing that causes the narrow GDR width in 120 Sn at low T. For this purpose we include the thermal 
pairing gap obtained from the modified-BCS theory [9T IToL ITU in the PDM 0. and carry out the calculations 

for the GDR width in 120 Sn at T < 5 MeV. 

The paper is organized as follows. Section I summarizes the main equations for the GDR including thermal pairing 
within the PDM and discusses in detail the physical assumptions of the PDM. Section II analyzes the results of 
calculations of GDR width, energy, and cross section for 120 Sn at finite temperature in comparison with the most 
recent experimental systematic. The paper is summarized in the last section, where conclusions are drawn. 
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II. MAIN EQUATIONS FOR HOT GDR WITHIN PHONON DAMPING MODEL 

A. Equations for GDR width and energy including thermal pairing 

The quasiparticle representation of the PDM including pairing has been already reported in Ref . ■ Therefore 
we discuss here only the final equations, which will be used in numerical calculations. According to this formalism 
the GDR width Tgdr is presented as the sum of quantal (Fq) and thermal (Ft) widths as 0, EL IS El 

Tgdr = r Q + T T , (la) 



T Q = 2-kFI Y^phfi 1 ~ n P~ nh)S(E GDn -Ep- E h ) , (lb) 

ph 

T T = 2nFi [v { J] 2 (ri s > ~ n s )6(E GTm - E s + E s .) , (lc) 

s>s' 

where (ss 1 ) = (pp 1 ) and (hh') with p and h denoting the orbital angular momenta j p and jh for particles and 
holes, respectively. The quantal and thermal widths come from the couplings of quasiparticle pairs [a p (g> a[]j,M and 
[aj ® ocs']lm to the GDR, respectively. At zero pairing they correspond to the couplings of ph pairs, [a P <8>a,h]LM, 
and pp (hh) pairs, [a J ® fl 8 ']z,Mi to the GDR, respectively (The tilde ~ denotes the time-reversal operation). The 
quasiparticle energies Ej — [(ej — ft,) 2 + A 2 ] 1 / 2 are found from the modified-BCS equations (39) and (40) of [Tl| . 
which determine the modified thermal gap A and chemical potential \x from the single particle energies Cj and particle 

number N. From them one defines the Bogoliubov coefficients Uj, Vj, and the combinations — u p Vh + v p Uh, and 



vL) = u.u. 



The quasiparticle occupation number rij is calculated as Q 



rij = — 

7T 



1 f°° n F (E)jj(E) <wE/T,-n-l f0 \ 

-I [E-E^M,(EW +1 2 (E) dE ' ME) = (e f +1) , (2) 

where Mj(E) is the mass operator, and jj(E) is the quasiparticle damping, which is determined as the imaginary part 
of the complex continuation of Mj(E) into the complex energy plane. These quantities appear due to coupling between 
quasiparticles and the GDR. Their explicit expressions are given by Eqs. (3) and (4) of Ref. respectively, in which 
E s is now the quasiparticle energy Ej . From Eq. J3J it is seen that the functional form for the occupation number rij is 
not given by the Fermi-Dirac distribution uf (Ej ) for non-interacting quasiparticles. It can be approximately to be so if 
the quasiparticle damping jj (E) is sufficiently small so that the Breit-Wigner-like kernel under the integration can be 
replaced with the <5-function. Equation also implies a zero value for rij in the grou nd state, i.e. rij(T = 0) = 0. In 
general, it is not the case because of ground-state corr elations (See. e.g Refs. [Sill They lead to n 3 (T = 0) ^ 
0, which should be found by solving self-consistently a set of nonlinear equations within the renormalized random- 
phase approximation (renormalized RPA). Within the RPA the equation for rij yields the approximate expression 
rij(T = 0) ~ X;jy/(2J+l)/(2j + l)[r^v 4) ] 2 , where is the RPA backward-going amplitude (J is the multipolarity). 

Since for collective high-lying excitations such as GDR one has < 1, we expect the value rij(T — 0) to be 

negligible. 

The GDR energy £"gdr is found as the solution of the equation 

uj - uj q - P(uj) = 0, (3) 
where uj q is the unperturbed phonon energy, and P(oj) is the polarization operator: 

P(E) = Fl V tf^ 1 -"*-"*) _ F| y ^?(ns-n s ,) 
y ' 1 V E-E p -E h 2 ^ E-E s +E s , W 

ph s>s 

Note that, in general, there are also backward-goi ng p rocesses leading to the terms ~ S(uj+E p +Eh) and S(uj+E s — E s >) 
as has been shown in Eqs. (14) and (15) of Ref. [l7|. However, as the maximum of these terms is located at negative 
energy uj = —(E p + Eh) < 0, and uj = —(E s — E s >) < 0, respectively, their contribution to the GDR, which is located 
at u = -Egdr S> 1 MeV, is negligible. Therefore these backward-going processes are omitted here. It is now easy to 
see that, at zero pairing A = 0, one has u p = 1, v p = 0, Uh = 0, Vh = 1 so that [w^] 2 = 1, [w^ s /'] 2 = 1. As for the 
single-particle occupation number fj, one obtains fh — 1 — Uh and f p — n p . The PDM equations for A = in Ref. 



smgle-particle occupation number jj, one obtains 
HQ are then easily recovered from Eqs. JU-JU 
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B. Assumptions of phonon damping model 

The PDM is based on the following assumptions: 

i) The matrix elements for the coupling of GDR to non-collective ph configurations, which causes the quantal width 
Tq l|lb|> . are all equal to F\. Those for the coupling of GDR to pp (hh), which causes the thermal width Tt (|Tc| . are 
all equal to F^. The assumption of a constant coupling strength is well justified when the width of a collective mode 
is much smaller than the energy range AE (of order of Ecrm) over which this mode is coupled to the background 
states (the so-called weak coupling limit discussed in Refs. [l3l 20}). 

ii) It is well established that the microscopic mechanism of the quantal (spreading) width Tq (jlbjl comes from 
quantal coupling of ph configurations to more complicated ones, such as 2p2h ones. The calculations performed in 
Refs. |2ll 22] within two independent microscopic models, where such couplings to 2p2h configurations were explicitly 
included, have shown that Tq depends weakly on T. The microscopic study in Ref. |20|. where an hierarchy of states 
of increasing complexity located around -Egdr is considered, has also confirmed the nearly constancy of Tq. It also 
indicated that the width of a collective vibration does not depend on the detailed coupling to the compound nuclear 
eigenstates. Therefore, in order to avoid complicate numerical calculations, which are not essential for the increase of 
Tgdr at T ^ 0, such microscopic mechanism is not included within PDM, assuming that Tq at T — is known. The 
model parameters are then chosen so that the calculated Tq and -Egdr reproduce the corresponding experimental 
values at T = (See below). 

Assumption (i) is satisfied for F\ since the quantal width Tq does not exceed 4.9 MeV. The PDM calculations in 
fact have shown that Tq decreases from 4.9 MeV at T = to around 2.5 MeV at T = 5 MeV for 120 Sn due to thermal 
effects in the factor 1 — n p — (See the dashed line in Fig. 1 (a) of for zero pairing). A similar trend was also 
observed in the microscopic calculations of Ref. |21 | . w here it was found that the GDR width at T — 3 MeV is in fact 
smaller that at T = (See Figs. 9 and 10 of Ref. [21| an the discussion therein). That's why Tq(T = 0) cannot be 
simply taken as a parameter uniformly added to what is calculated for ^ at T ^ since a width Tq(T) = Tq(T = 0) 
would lead to a larger value for the total width Tgdr <|lafl at higher T, worsening the agreement with the data. 

Assumption (i) becomes poor for F2 at T > 3 MeV, when the thermal width Tt is larger than 10 MeV (See the 
dotted line in Fig. 1 (a) of Ref. 0). Within such a large width one expects a considerable change of the level density 
of background states. To be quantitatively precise, one needs to use a self-consistent theory for the strength function, 
which includes the coupling to doorway states as in Ref. |2jj. Such theory is valid for any ratio of Fgdr/AE. 
However, from assumption (ii) it also follows that the increase of Tgdr is now driven mainly by the thermal width Tt 
due to the factor n s i —n s . The change of rij implies a change of the quasiparticle entropy S qp . The latter is ultimately 
related to the change of the level density of background states within the realistic mean-field basis. This can be seen 
as follows. The complexity of the background states is measured by the information entropy Sj n f of individual wave 
functions, which reflects the complicated relationship between the eigenbasis and representation basis. Meanwhile, 
the thermodynamic entropy Sth of the total system is directly determined by its statistical weight Q(E) — p(E)8(E) 
as Sth = hif2(E), where p(E) is the level density. In a situation with incomplete information, such as in the statistical 
description of hot nuclei considered here, individual compound systems are replaced with a grand canonical ensemble 
of nuclei in thermal equilibrium. The probability for a quantum system to have a given eigenenergy is determined by 
the density matrix T> rather than by a pure wave function. The expectation value (O) of an observable O is given 
as the statistical average over the grand canonical ensemble (O) = Tr(PO), which is derived from the maximum of 
the thermodynamic entropy Sth = — Tr(2?lnD). The modern shell- model calculations in Ref. |13| have shown that 
these three apparently different entropies, S qp , S; n t, and Sth, behave very similarly for the majority of states in the 
realistic mean field consistent with residual interactions (See column II of Fig. 56 in Ref. [Ij] and the discussion 
therein). Significant differences between them take place only when the residual interaction beyond the mean field is 
very weak (See column I of Fig. 56 in |13l) or when the quasiparticle mean- field is absent (Sj n f reaches its chaotic 
limit) (See column III of Fig. 56 in |l3j|). This might be the reason why the numerical results performed so far 
within the zero-pairing PDM 0, 0, using assumptions (i) and (ii) fit fairy well the experimental systematic for the 
GDR width including the existing data on the width saturation at T > 3 MeV. This is partly also due to the large 
experimental errorbars for the extracted GDR width at high T (See, e.g. Refs. (2^, 0])- Therefore, at T > 3 MeV 
the results obtained under these assumptions should be considered as qualitative. This does not affect our present 
study of the pairing effect as the latter is significant only at low temperature (T < 2.5 MeV). 

Within assumptions (i) and (ii) the model has only three T-independent parameters, which are the unperturbed 
phonon energy uo qi F\, and F2. The parameters ui q and F\ are chosen so that after the p/i-GDR coupling is switched 
on, the calculated GDR energy -Egdr and width Tgdr reproduce the corresponding experimental values for GDR 
on the ground-state. At T ^ 0, the coupling to pp and hh configurations is activated. The F2 parameter is then 
fixed at T = so that the GDR energy -Egdr does not change appreciably with varying T. The values of the PDM 
parameters for 120 Sn are given in |4j for the zero-pairing case. 

In Ref. Q we have presented an argument that, in our opinion, the effect due to coupling of GDR to noncollective 
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ph, pp and hh configurations at T 7^ within the PDM is tantamount to that of thermal shape fluctuations. This 
has been demonstrated by expanding the coupling of GDR phonon to noncollective ph, pp, and hh configurations into 
couplings to different multipole fields (See pages 437 and 438 of Ref. |4|). In this expansion the pp (hh)-p&ir operator 
B\ s , = a\a s > is expanded in terms of the tensor products of two ph-p&ir operators (See mappings (22) and (23) of Ref. 
|19|'). Each p/i-pair operator B^ h = a^a^ can be then expressed in terms of the RPA phonon operators Q q and Q q 
with RPA amplitudes X q ph and Y* h . This leads to Eq. (2.43) in Ref. for the part of the PDM Hamiltonian, which 
describes coupling between the phonon and single-particle fields. Hence if and Q q are GDR phonon operators, 
{Q^iQqi} an d {Q\ 2 jQq 2 } m Eq. (2.43)ofRef. Q can have the moment and parity equal to (1~, 2+) , (2 + ,3~),etc. to 
preserve the total momentum A 71 " = 1~. Therefore, although F^} are dipole matrix elements, the amplitudes X^ h and 
(i — 1,2) can be calculated microscopically, using the dipole-dipole, quadrupole-quadrupole, octupole-octupole, 
etc. components of residual interaction. This means that coupling to pp and hh configurations already includes in 
principle the coupling to different multipole- multipole fields via multiphonon configuration mixing at T ^ 1 . 

cascade calculations using PDM strength functions have produced the GDR shapes in good agreement with the 
experimental data for 120 Sn at T > 1.54 MeV (See Fig. 2 of Ref. 7]). The discrepancy at T < 1.54 MeV is due 
to omission or improperly inclusion of pairing at low T, which is now studied in the present work. The splitting of 
GDR into two peaks is also clearly seen in the PDM calculations for 6 Sn These evidences are in favor of the 
argument discussed above and in Ref. Q. 

Nonetheless, we recognize that the issue of whether coupling to pp and hh configurations at T ^ is a microscopic 
(although indirect) interpretation of thermal shape fluctuations within the PDM is still not settled. The question on 
whether thermal shape fluctuations need to be included additionally within PDM or not remains to be investigated. 
The thermal shape-fluctuation model calculates the time-correlation function of the GDR by replacing the micro- 
canonical ensemble with an ensemble of macroscopic variables, which are the deformation parameters (/3, 7) in the 
body-fixed (principal axes) frame of reference (intrinsic frame). It parametrizes a priori the dipole correlation tensor 
by a frequency — Eq exp [— yS/ 47T/3 cos(7 + §7rfc)] and a width equal to = Tq^E^/ Eq) 1 - 8 along each ith semiaxis 
[2$ . The effect of thermal fluctuations in this model is included via fluctuations of shapes by employing the macroscopic 
Landau theory of phase transitions |27j. Therefore, an explicit inclusion of thermal shape fluctuations in the PDM 
will bring in additional degrees of freedom, which increase the number of parameters of the model. At the same 
time, in a way similar to that mentioned in the footnote above, this will certainly require a renormalization of the 
existing parameters of the PDM to restore the agreement with the experimental systematic. While this issue is left 
open for future study, it does not affect the study of the role of thermal pairing in the present paper, since, as will be 
seen below, in order to describe the GDR width at low T, thermal pairing is necessary to be included in any model, 
whether it is the PDM or thermal shape-fluctuation one. 



III. ANALYSIS OF NUMERICAL CALCULATIONS 

A. Role of thermal pairing gap at low temperature 

Shown in Fig. ^is the T dependence of the neutron pairing gap A„ for 120 Sn, which is obtained from the modifled- 
BCS equations [Uli E3 

using the single-particle energies determined within the Woods-Saxon potential at T = 0. 
They span a space from ~ -40 MeV up to ~ 17 MeV including 7 major shells and lj'15/2, l*n/2i an d 1^17/2 levels. The 
pairing parameter G v is chosen to be equal to 0.13 MeV, which yields A(T = 0) = A(0) ~ 1.4 MeV. In difference with 
the BCS gap (dotted line), which collapses at T c ~ 0.79 MeV, the gap A (solid line) does not vanish, but decreases 
monotonously with increasing T at T > 1 MeV resulting in a long tail up to T ~ 5 MeV. This behavior is caused by 
the thermal fluctuation of quasiparticle number, 5M 2 = JV <5A/J, where SMJ — nj( \ —rij) is the quasiparticle- number 
fluctuation on j-th orbital. The latter is incorporated as SAfj in the modificd-BCS gap A (Sec the last term at the 
r.h.s of Eq. (39) of Ref. [HI). 



1 In Ref. 1251 a version of PDM, which explicitly includes coupli ng t o two-phonon configurations in the second order of the interaction 
vertex, was proposed. The expressions (2.17) and (2.18) of Ref. 25| derived for the polarization operator show that the width increase 
is still driven mainly by the factor (n s — n s i). The calculations with the GDR coupled to the first quadrupole phonon 2^~ required the 

energy and the ratio r = F\ /F^ (i = 1,2) to be fixed as additional parameters. A similar quality for the description of the 
experimental data for the temperature dependence of the GDR width has been restored after a reduction of the dipole matrix elements 
= F 1 and F 2 (1) = F 2 . 
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To analyze the qualitative effect of thermal pairing on the GDR width we plot the usual single-particle occupation 
number fj — ujrij + u?(l — fij) with rij obtained within the modified-BCS theory ( r )j(E) — 0) as a function of single- 
particle energy ej for the neutron levels around the chemical potential in Fig. [21 It is seen that, in general, the pairing 
effect always goes counter the temperature effect on fj, causing a steeper dependence of fj on Ej. Decreasing with 
increasing T, this difference becomes small at T > 3 MeV. Since a smoother fj enhances the pp and hh transitions 
leading to the thermal width Tt , pairing should reduce the GDR width, and this reduction is expected to be stronger 
at a lower T, provided the GDR energy -Egdr is the same. A deviation from this general rule is seen at very low T ~ 
0.1 MeV, where the temperature effect is still so weak that fj obtained at A ^ (solid line) is smoother than that 
obtained at zero pairing (dotted line). 

To get an insight into the detail of the change of GDR width at low T we show in Fig. [21 tlie combinations [u^] 2 

and [VtL] 2 of the Bogoliubov coefficients u p , Vh, and iv together with the factors (1 — n p — rih) and (rih — rih>) as 
well as their zero-pairing counterparts, fh — f P and fh> — fh for the particle p = 2j 7 / 2 , hole h = 2d 3 / 2 , and b! = ld 5 / 2 
orbits as functions of T. The hole orbit h = 2d$/2 is located just below the chemical potential. Therefore the pairing 

effect is strongest for the ph and hh configurations involving this orbit. This figure shows a sharp increase of [u^'] 2 
and [v^] 2 at T < 1 - 2 MeV due to a steep slope of the pairing gap, showing a strong pairing effect. At very low T, 
the numerator of the polarization operator (0J is close to [u^] 2 (E p + Eh) because 1 — n p — rih — 1, while the thermal 
part ~ (n s — n s >) ~ 0. This value is equal to 3.67 MeV at T = 0.1 MeV, which is smaller than (eh — e p ) = 4.63 MeV 
obtained at A = 0. The denominator of (0} is also smaller than that obtained at A — because E p + Eh > eh — £p 
due to the gap. Therefore, at very low T, pairing may lead even to a smaller GDR energy. On the other hand, as 
[ u pt^ 2 Hi — n p — rih) and [v[~}} 2 (n s — n s >) are also smaller than fh — f p and f s > — f s , respectively, the competition of 
these effects in Eq. (1) can result in a. larger width in the very low-T region. As T increases, the factor 1 — n p — •;»./, 

decreases while [w^] 2 increases to reach 1 at T > 2 MeV because of the decreasing gap. This leads to the decrease of 
the quantal width Tq . At the same time, coupling to pp and hh configurations starts to contribute due to the factor 

(n s — n s i). The combination [«£,^] 2 a l so increases with T and reach 1 at T > 3 MeV. As the result, thermal width 
Tt starts to give an increasing contribution with T. However, as compared to their zero-pairing counterparts, fh — f p 
and f s i — f s , the decrease of the quantal part ~ (1 — n p — rih) and increase of the thermal part ~ (n s — n s <) are much 
more moderate with increasing T up to 1 MeV. On the contrary, at 1 < T < 3 MeV the decrease (1 — n p — rih) and 
increase of (n s — n s >) with increasing T are steeper than their counterparts at A = 0. At T > 3 - 4 MeV the total 
width approaches the saturation because of the dominating contribution of Tt, which ceases to increase due to the T 
dependence of n s — n s < shown in Fig. [31 (b) 



B. Temperature dependence of GDR width and energy 

The GDR width Tgdr and energy -Egdr for 120 Sn were calculated from Eqs. (JTJ and ©, respectively, using the 
same set of PDM parameters iv q , Fi, and F2, which have been chosen for the zero-pairing case 03 (set A). The 
effect of quasiparticle damping is included in the calculations by using Eq. (J2J. The results are shown as the thin 
solid lines in Fig. 01 As seen from Fig. |U(b), the oscillation of GDR energy -Egdr with varying T occurs within 
the range of ~ ±1.5 MeV, which is wider compared with that obtained neglecting pairing. The latter is almost 
independent on T (dashed line in Fig. 0] (b)). As expected from the discussion above, the GDR energy £gdr at 
T = 0.1 MeV drops to 14 MeV, i.e. by 1.4 MeV lower than the GDR energy measured on the ground state. The 
GDR width increases to 5.3 MeV compared to 4.9 MeV on the ground state as shown in Fig. 01 (a). At < T < 
0.5 MeV, the above-mentioned competition between the decreasing quantal and increasing thermal widths makes the 
total width decrease first to reach a minimum of 3.4 MeV at T ~ 0.2 MeV then increase again with T. At T > T c 
the width only increases with T. At 1 < T < 3 MeV the GDR width obtained including pairing is smaller than the 
one obtained neglecting pairing (dashed line in Fig. El (a)), but this difference decreases with increasing T so that at 
T > 3 MeV, when the gap A becomes small, both values nearly coincide. This improves significantly the agreement 
with the experimental systematic at 1 < T < 2.5 MeV. In order to have the same value of 4.9 MeV for the GDR 
width at T ~ 0, we also carried out the calculation using slightly readjusted values F[ — 0.96Fi and F' 2 = 1.03^2 
while keeping the same u> q (set B). The result obtained is shown in the same Fig. 01 as the thick solid lines. The GDR 
energy -Egdr moves up to 16.6 MeV at T = 0.1 MeV and at T c < T < 1.2 MeV in agreement with the value of 16.5 ± 
0.7 MeV extracted at T — 1 MeV in Ref. The width at T = 1 MeV also becomes slightly smaller, which agrees 
quite well with the latest experimental point [l6|. At T > 1 MeV the results obtained using two parameter sets, 
A and B, are nearly the same. The effect of quantal fluctuations SN 2 of particle number within the BCS theory at 
T = 0, however, is neglected in these results. To be precise, this effect should be included using the particle-number 
projection method at finite T. However, the latter is so computationally intensive that the calculations were carried 
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out so far only within schematic models (See, e.g. |3j|), or one major shell for nuclei with A < 60 as in the Shell- 
Model Monte-Carlo method . Therefore, for the limited purpose of the present study, assuming that 8N 2 3> 1, 
we applied the approximated projection at T — proposed in Ref. |33| . which leads to the renormalization of the 
gap as A(T) = [1 + 1/SN 2 ]A with SN 2 = A(0) 2 YV(j + V 2 )/[(ej - Z 2 ) 2 + A(0) 2 ] 0- This yields A(T = 0) ~ 1.5 
MeV (5N ~ ± 4). The PDM results obtained using the gap A and the parameter set B are shown in the same Fig. 
0]as the thick dotted lines. The GDR width becomes 5 MeV with -Eqdr = 15.3 MeV at T — in good agreement 
with the GDR parameters extracted on the ground state. It is seen that the fluctuation of the width at T < 0.5 MeV 
is largely suppressed by using this renormalized gap A. For comparison, the predictions by two versions of thermal 
shape-fluctuation model 0, Eg are also plotted in Fig. 0] (a) as the dash-dotted and thin dotted |2^| lines. It 
is seen that these predictions, in particular the one given by the phenomenological version in Ref. |28j| . significantly 
overestimate the GDR width at low temperature T < 1.3 MeV. The predicted overall increase of the width is not 
as steep as the experimental systematic and the PDM prediction. The curvature of the trend is also opposite to 
the experimental one and that given by the PDM. It is, therefore, highly desirable to see how the prediction by the 
thermal shape-fluctuation model would change by taking into account the effect of thermal pairing gap discussed in 
the present work in combination with the use of a specific Hamiltonian to calculate every quantity 35] . 



C. Effect of thermal pairing on GDR cross section at low temperature 



Shown in Fig. [3]are GDR cross-sections obtained for 120 Sn using Eq. (1) of Ref. 0. The experimental cross-section 
are taken from Fig. 2 of Ref. 0. They have been generated by cascade at excitation energy E* — 30 and 50 MeV, 
which correspond to r max = 1-24, and 1.54 MeV, respectively. The theoretical cross-sections have been obtained using 
the PDM strength function Sgbr(E 7 ) from Eq. (2.22) of Ref. [f| at T — T max . This is the low temperatures region, 
at which discrepancies are most pronounced between theory and experiment. (A divided spectrum free from detector 
response at T = 1 MeV is not available in Ref. [3). From this fi gure it is seen that thermal pairing clearly offers a 
better fit to the experimental line shape of the GDR at low temperature. As has been discussed in Ref. '7\, for an 
absolute comparison, the PDM strength functions for all decay steps starting from T max down to T ~ MeV should 
be included in the cascade to generate a divided spectra, which can be directly compared with the experimental 
ones. It is our wish that such calculations can be carried out in collarboration with the authors of fl6| in the near 
future. 



IV. CONCLUSIONS 



In this paper we have included the pairing gap, determined within the modified-BCS theory 0,^3, in the PDM 
to calculate the width of GDR in 120 Sn at T < 5 MeV. In difference with the gap given within the conventional BCS 
theory, which collapses at T c ~ 0.79 MeV, the modified-BCS gap never vanishes, but monotonously decreases with 
increasing T up to T ~ 5 MeV. The results obtained show that thermal pairing indeed plays an important role in 
lowering the width at T < 2 MeV as compared to the value obtained without pairing. This improves significantly 
the overall agreement between theory and experiment, including the width at T = 1 MeV extracted in the latest 
experiment |lq|. 



Acknowledgments 

The numerical calculations were carried out using the FORTRAN IMSL Library 3.0 by Visual Numerics on the 
Alpha server 800 5/500 at the Division of Computer and Information of RIKEN. 



[1] M.N. Harakeh and A. van der Woude, Giant resonances - Fundamental high-frequency modes of nuclear excitation (Claren- 
don Press, Oxford, 2001). 

[2] W.E. Ormand, P.F. Bortignon, and RA. Broglia, Phys. Rev. Lett. 77, 607 (1966); W.E. Ormand et al, Nucl. Phys. A 614, 
217 (1997). 

[3] N.D. Dang and A. Arima, Phys. Rev. Lett. 80, 4145 (1998). 

[4] N.D. Dang and A. Arima, Nucl. Phys. A 636, 427 (1998) 

[5] N. Dinh Dang, K. Tanabe, and A. Arima, Nucl. Phys. A 645, 536 (1998). 



7 



[6] G. Gervais, M. Thoennessen, and W.E. Ormand, Phys. Rev. C 58, R1377 (1998). 

[7] N. Dinh Dang, K. Eisenman, J. Seitz, and M. Thoennessen, Phys. Rev. C 61, 027302 (2000). 

[8] L.G. Moretto, Phys. Lett. B 40, 1 (1972). 

[9] N. Dinh Dang and A. Arima, Phys. Rec. C 68, 014318 (2003). 

[10] N. Dinh Dang and V. Zelevinsky, Phys. Rev. C 64, 064319 (2001). 

[11] N. Dinh Dang and A. Arima, Phys. Rev. C 67, 014304 (2003). 

[12] N.D. Dang, P. Ring, and R. Rossignoli, Phys. Rev. C 47, 606 (1993). 

[13] V. Zelevinsky, B.A. Brown, N. Frazier, and M. Horoi, Phys. Rep. 276, 85 (1996). 

[14] A. Volya, B.A. Brown, and V. Zelevinsky, Phys. Lett. B 509, 37 (2001). 

[15] N. Dinh Dang, K. Tanabe, and A. Arima, Nucl. Phys. A 675, 531 (2000). 

[16] P. Heckman et al., Phys. Lett. B 555, 43 (2003). 

[17] N.D. Dang, V.K. Au, T. Suzuki, and A. Arima, Phys. Rev. C 63, 044302 (2001). 

[18] K. Hara, Prog. Theor. Phys. 32, 88 (1964); K. Ikeda, T. Udagawa, and H. Yamamura, Prog. Theor. Phys. 33, 22 (1965); 

D.J. Rowe, Phys. Rev. 175, 175 (1968); J. Dukelsky and P. Schuck, Nucl. Phys. A 512, 466 (1990). 
[19] F. Catara, N. Dinh Dang, and M. Sambataro, Nucl. Phys. A 579, 1 (1994). 

[20] B. Lauritzen, P.F. Bortignon, R.A. Broglia, and V.G. Zelevinsky, Phys. Rev. Lett. 74 (1995) 5190. 

[21] P.F. Bortignon, R.A. Broglia, G.F. Bertsch, and J. Pacheco, Nucl. Phys. A 460, 149 (1986). 

[22] N. Dinh Dang, Nucl. Phys. A 504, 143 (1989). 

[23] A. Bracco et al, Phys. Rev. Lett. 62, 2080 (1989). 

[24] J.J. Gaardhoje et al, Phys. Rev. Lett. 59, 1409 (1987) 

[25] N. Dinh Dang, K. Tanabe, and A. Arima, Phys. Rev. C 58, 3374 (1998). 

[26] N. Dinh Dang, A. Ansari, and A. Arima, preprint RIKEN-AF-NP-361 (2000), N. Dinh Dang, Nucl. Phys. A 687, 253c 
(2001). 

[27] L.D. Landau and E.M. Lifshitz, Course of Theoretical Physics, Vol. 5: Statistical Physics (Moscow, Nauka, 1964). 

[28] D. Kusnezov, Y. Alhassid, and K.A. Snover, Phys. Rev. Lett. 81, 542 (1998) and references therein. 

[29] M.P. Kelly et al., Phys. Rev. Lett. 82, 3404 (1999). 

[30] B.L. Berman, At. Data Nucl. Data Tables 15, 319 (1975). 

[31] R. Rossignoli, N. Canosa, and P. Ring, Phys. Rev. Lett. 80, 1853 (1998). 

[32] Y. Alhassid, S. Liu, and H. Nakada, Phys. Rev. Lett. 83, 4265 (1999). 

[33] I.N. Mikhailov, Sov. Phys. JETP 18, 761 (1964). 

[34] N.D. Dang, Z. Phys. A 335, 253 (1990). 

[35] A. Ansari, N. Dinh Dang, and A. Arima, Phys. Rev. C 62, 011302(R) (2000), Phys. Rev. C 63, 024310 (2000). 



8 



figl.gif 

FIG. 1: Neutron pairing gap as a function of T. Solid and dotted lines show the modified-BCS gap A and BCS gap, respectively. 

fig2.gif 

FIG. 2: Single-particle occupation number fj as a function of ej for the neutron levels around the chemical potential at T = 
0.1, 1, and 3 MeV. Results obtained including and without pairing are shown by solid and dotted lines, respectively (A thicker 
line corresponds to a higher T). The horizontal dashed line at ~ -6 MeV shows the chemical potential at A = and T = 0. 



fig4.gif 

FIG. 4: GDR width Tgdr (a) and energy Sgdr (b) as functions of T for 120 Sn. The black solid lines show the PDM results 
obtained neglecting pairing (In (a) it is the same as the solid line with diamonds from Fig. 1 (a) of Ref. 3j). The green and red 
lines are the PDM results including the gap A, which are obtained using the parameter sets A and B, respectively. The blue 
lines are the PDM results including the renormaliz ed g ap A (See text). Solid circles are the low-T data from Ref. ("ill . Crosses 
and open triangles in (a) are from Fig. 4 of Ref. |2S|. The corresponding GDR energies decrease from 16 MeV to 14.5 MeV 
with increasing T as shown in the shaded rectangle in (b). Solid upside-down triangles are data from Ref. [2Sll . Open squares 
and stars are high-T data for 110 Sn from Refs. |23| and |24|, respectively. Data at T = are for GDR built on the ground state 
of tin isotopes with masses A = 116 - 124 from Ref. [3C§- The predictions by two versions of the thermal shape-fluctuation 
model are shown in (a) as the dash-dotted and thin dotted |2S| lines, respectively. 



fig5.gif 

FIG. 5: Experimental (shaded areas) and theoretical divided spectra obtained without pairing (black solid lines) and including 
the gap A (red lines) as for the red lines in Fig. 2] 



fig3.gif 

FIG. 3: Combinations of Bogoliubov coefficients (a) and quasiparticle occupations numbers (b). In (a): The thin and thick 
lines show [w^'] 2 and [i>^~/] 2 , respectively, for the orbits p — 2j 7 / 2 , h = 2d 3 / 2 , and hi = lds/2- I n (b)i The corresponding 
factors 1 — rip — rih (thin line) and nh — n h i (thick line) are shown. The factors fh — f P and f h i — fh for the zero-pairing case 
are shown as the dotted and dashed lines, respectively. 
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